The Latent Heat of Single Flavor Color Superconductivity in a Magnetic Field 
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We calculate the energy release associated with first-order phase transition between different types 
of single flavor color superconductivity in a magnetic field. 
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At sufficiently high baryon density, a hadronic mat- 
ter will be squeezed into a degenerate system of quarks, 
which becomes a color superconductor(CSC) at suffi- 
ciently low temperature [lj]. 

The most plausible place in the universe where the 
color superconductivity may be implemented is the inte- 
rior of compact stars. They have masses about 0.5A/q — 
2.7Mq and are believed to have radii of order 10km. 
Near the surface their density is around normal nu- 
clear density and raises to several times higher in the 
core region. There a quark matter of chemical potential 
^ ~ 400 - 500MeT^ may form. 

The color superconductivity can influence the observ- 
able signature of a compact star in several ways. Firstly, 
it will impact on the quark matter equation of state. Sec- 
ondly, the energy gap, as well as the Goldstone and/or 
Higgs modes associated to the long range interaction the 
transport properties of the star. Thirdly, the presence of 
strong magnetic field renders the superconducting phase 
transition first order. The latent heat release at the tran- 
sition may lead to observable energy burst which will be 
discussed in this report. 

While the color-fiavor-locked(CFL) phase, which pairs 
different quark flavors, is the ground state at ultra high 
baryon density (chiral limit), the situation is far more 
complicated at the moderate density in a compact star. 
The mass of s quarks and the electrical neutrality con- 
straint create a substantial mismatch of the Fermi mo- 
menta among different flavors which in turn reduces the 
available phase space of Cooper pairing. The single-flavor 
CSC (pairing within each flavor), however, is free from 
such a limitation. The smallness of the single-flavor gaps 
because of the nodal lines of the pairing force, makes the 
single-flavor phases relevant for late age stars, when the 
temperature is sufficiently low. Among the four canonical 
phases, the spherical color-spin-lock(CSL) and nonspher- 
ical planar, polar and A0-(1], CSL is the most favored one 
in the absence of a magnetic field. 

The presence of a magnetic field offsets the energy 
balance among the four canonical single flavor pairings. 
Only the spherical CSL phase has electromagnetic Meiss- 
ner effect [2j and nonspherical phases: polar, A and pla- 
nar phases only shield part of the magnetic field. So 
if a quark matter cools down through the critical tem- 
perature of the single flavor paring in a magnetic filed, 
forming CSL state will cost extra work to exclude mag- 
netic fluxes from the bulk. Therefore, the magnetic con- 



tribution to the free energy may favor the non-spherical 
states. In a previous work[5j, we have explored the con- 
sequences of the absence of the electromagnetic Meissner 
effect in a nonspherical CSC phase of single-flavor pair- 
ing and have obtained the phase diagram with respect 
to the temperature and the magnetic field. Computing 
the latent heat released acrossing the phase boundaries 
is the main subject of the present work. We shall employ 
the same approximation inQ by ignoring the masses of 
pairing quarks. 

We work with a NJL-likc effective action whose La- 
grangian density reads jg|: 

C = 4>{-^ v d v + ^)vb + Gvb lv T l ^ lv T 1 ^ (1) 

with T l = iA , G > an effective coupling and introduce 
the condensate: 



(2) 



where tp is the quark field, %j)c = 72"0* is its charge con- 
jugate, A c with c = 2, 5, 7 are antisymmetric Gell-Mann 
matrices and r c is a 4 x 4 spinor matrix. We may choose 
T 5 = T 7 = for the polar and A phases, T 2 = for the 
planar phase but none of r c 's vanishes for CSL phase. 
We find the pressure of each flavor under mean field ap- 
proximation: 

P= - J^- A -£ k )-^E( lt -' 4 -^-Ml) 
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where E k = y/Qe - n) 2 + A 2 f 2 (8) with 8 the angle be- 
tween the momentum k and a prefixed spatial direction 
and A the gap parameter. The function f{9) is given by 



1, for CSL phase 



f(0) = 



— (1 + cos 2 8), for planar phase 
sin 8, for polar phase 
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cos -. for A phase 
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Maximizing the pressure with respect to A, we obtain the 
gap equation — 0, which determines the temper- 

ature dependence of the gap, A(T). Substituting A(T) 

back to ©, WC find that P n < P A < Ppolar < -Pplanar < 

■PcsL up to the transition temperature T c . 

The diquark condensate ([2]) for CSL breaks the gauge 
symmetry SU(3) C x U(l) em completely. A nonspherical 
condensate, however, breaks the gauge symmetry par- 
tially and the Meissner effect is incomplete. Among the 
residual gauge group, which is the gauge transformation 
that leaves the diquark operator inside @ invariant, there 
exists a U(l) transformation, ip — > e~' A80_49 *V with q 
the electric charge of tp, 6 = —2y/3q(f> for the polar and 
A phases and 6 = 4y/Sq<p for the planar phase. The 
corresponding gauge field, A^ is identified with the elec- 
tromagnetic field in the condensate and is related to the 
electromagnetic field A and the 8-th component of the 
color field A 8 in the normal phase through a U(l) rota- 
tion 

Aft = Apcosj- ApSW.'y 

V M = A M sin 7 + A® cos 7 (5) 

where the mixing angle 7 is given by tan7 po i arj A = 
2^/3q(e/g) and tan7 P i anar = 4^/3q(e/g) with g the QCD 
running coupling constant. The 2nd component of (J5]) 
V = because of the Meissner effect and thereby imposes 
a constraint inside a nonspherical CSC, A 8 = —A^ tan 7, 
which implies that: 

B 8 = -Btan 7 (6) 

The thermal equilibrium in a magnetic field Hz is de- 
termined by minimizing the Gibbs free energy density, 

G = -P+ l -B 2 + l -Y.(B l f-BH (7) 
1=1 

with respect to A, B and B l . Ignoring the induced mag- 
netization of quarks, the pressure P is given by with 
A given by the solution of the gap equation. For a non- 
spherical CSC pairing, the minimization with respect to 
B and B l is subject to the constraint ©. For a hypo- 
thetical quark matter of one flavor only, we find that 

Smin. j = -Pi - \vjH 2 (8) 

with j = n, CSL, polar, A and planar, where rj n = 1, 
?7csl = and rjj = cos 2 7, for a nonspherical CSC. The 
phase corresponding to minimum among £/ m i n 's above 
wins the competition and transition from one phase to 
another is first order below T c . 

For a quark matter of different flavors with Cooper 
pairing within each flavor, the Gibbs free energy density 
can be written as 

g = -P-h n H 2 . (9) 



where P is the total pressure of all flavors. We have 
r] = 1 if all flavors are normal and r\ = if all flavors 
are CSL. If one flavor is in a nonsherical state and others 
are normal, 77 = cos 2 7 with 7 the mixing angle of the 
nonspherical phase. If more than one flavors are in non- 
spherical phases, r\ — cos 2 7 with 7 their common mixing 
angle. If their mixing angle were different, we would end 
up with B = £> 8 = 0, in order to compromise the con- 
straints ([6]) of all nonspherical states involved, making 
them less favored than CSL. It was shown in Q that 
only four phases, I-IV(shown in Table. I and Fig.l), need 
to be considered in both two and three flavor quark mat- 
ters. There is a first order phase transition from one of 
them to another for < T < T c . 
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FIG. 1: The H-T phase diagram for two flavors and three 
flavors, where the reference magnetic field Ho = /j.Ao/tt with 
A the CSL gap at T = 0. 

The border between two phases are determined by the 
equation: 

H 2 H 2 

Pa+Va— = PP+VP— (10) 

with the subscripts a and f3 labeling the four phases I-IV 
and the density of the latent heat released from the phase 
a to the phase j3 reads 

Q a p=T[S a (T)-S p (T)], (11) 

The entropy density S(T) is given by 

(dp\ _ fdp\ fdp\ fdA\ _ /ap\ 
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TABLE I: The possible phases involving spin-one CSC of a quark matter with two flavors or three flavors in a magnetic field. 





I 


ii 


III 


IV 


2 flavor 


CSL U , CSLd 


(polar) u , (planar) d 


(normal) u , (polar)d 


(normal) u , (normal)d 


3 flavor 


CSL U , CSLd.s 


(polar) u , (planar) d , s 


(normal) u , (polar) d , s 


(normal) u , (normal)d, s 



where the gap equation is employed in the last step. 

The temperature dependence of the thermodynamic 
quantities, P and S can be written in a parametric form. 
In terms of the parameter t — , the gap equation 

h(t) with 



takes the form In 



h(t) = [ d9 sm6 f 2 {6) 
Jo 

io°° dx 



(12) 



^/x 2 +t 2 f 2 (6)[e^ x2+t2f2 ^ + 1] 



It follows that 
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(13) 



Introducing P s - P n = p(t)^r and S s - S n = 
m 2t ^° f df/^jf^ with s lableling different pairing states 
and A = Aqsl(O), we have 



p(t) = e- 2h ^ [1 + 2h(t) + 4^ - (14) 



,ffW_27r_ 2 
t 2 3t 2 



with 

g{t)= ( ddsmO f dx\n[e-^ x2+t2f2{ ^ + 1] (15) 
Jo Jo 

and the curves P(T) and S(T) may be plotted paramet- 
rically without solving the gap equation for T > 0. As 
we did in 

Numerically, we identify Ao with that of the one-gluon 
exchange 0, 0] with pL = 5Q0MeV and a s = f^r = 1 as 
a calibration of the parameters in (fT]), We find Ao — 
0.238MeV, H ~ 5.44 x 10 14 G for two flavors, and 
A ~ 0.0864MeU, H ~ 1.97 x 10 14 G for three flavors. 
As is shown in Fig.l, the nonspherical phases occupy a 
significant portion of the H-T phase diagram for a mag- 
nitude of the magnetic field of order 10 15 G. This magni- 
tude of the magnetic field is plausible in a compact star. 

The temperature dependence of the entropy density 
differences between all single flavor CSC phases and that 
of the normal phase under the same chemical potential 
fi = 500AfeV is shown in Fig. 2. The differences vanish 
at T — since all entropies are zero. They also van- 
ish at T — T c since the transition is second order there. 
The entropy density of a CSC phase is lower than that 
of the normal phase because of the long range order, so 
the differences are all negative for < T < T c . The 
latent heat densities across the phase boundaries in the 
Fig.l are displayed in the Fig. 3, where adjustment of the 
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FIG. 2: The entropy density differences between single flavor 
CSC phases and normal phase as a function of the tempera- 
ture. The unit is erg/ K/km i . 



chemical potential of each flavor is made to fulfill the re- 
quirement of the charge neutrality. The latent heat gets 
a peak value in about 0.67T C -0.72T C for all curves. 

In the natural unit we are using, the entropy is dimen- 
sionless and the unit of latent heat can be MeV , the unit 
of latent heat density is MeV 4 = 2.088 x 10 41 erg/km 3 . 
The units have been transformed to International Sys- 
tem of Units in Fig. 2 and Fig. 3. For the biggest latent 
heat density we calculated in three flavor phase, from 
IV to III, this value is about 6.04 x 10 44 erg/km 3 . Then 
for a compact star with the radii of quark matter equal 
to R(km), the energy release is about 2.53 x 10 45 R 3 erg. 
Since R < 10km typically, the latent heats released as 
the star cools through the phase boundaries we showed 
are smaller than the typical energy release of 7 ray burst 
of order 10 51 to 10 54 erg[8|, Q. But they may contribute 
to weaker energy bursts such as the X-ray radiation at 
the later stage of a compact star[Toj. 

Finally, we would like to comment on the approxima- 
tions employed. The two flavor case is certainly unreal- 
istic since it assumes a large mass of s quarks, m s >> p. 
For three flavors, our approximation requires that m s to 
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FIG. 3: The latent heat density dependence of tempera- 
ture. We know that S s = S n at T c and Q Q/3 = T[S a (T) - 
Sp(T)], so latent heat also vanish at To and T c . The unit of 
Q is erg /km 3 . 



be sufficiently smaller than [i. On the other hand m s can- 
not be too small in order for the single flavor CSC to com- 
pete with multi-flavor pairings under the Fermi momen- 
tum mismatch, which requires m s to be of the same or- 
der of magnitude as \J fiA(2SC) or greater, with A(25C) 
the 2SC gap at zero temperature without mismatch. The 
phase structure of 2SC in a magnetic field has been stud- 
ied in|ll|. Since A(2SC) is about 5-10 times lower than 



/x for the quark matter we considered, the approximation, 
though marginal, may not be too crude as is suggested by 
our recent analysis with a realistic m s , 150MeV, in the 
single flavor pairing dynamics [l^. For the highest mag- 
netic field in the phase diagram of Fig.l, H ~ 5 x 10 15 G, 
which implies e_ff/(/i 2 ) ~ 6 x 10~ 5 . Therefore it is legit- 
imate to ignore the impact of the magnetic field on the 
pairing dynamics, unlike the situation considered in (ljjj . 

Throughout the paper we have assumed that the whole 
volume of the quark matter core of a compact star un- 
dergoes a phase transitions at the same time. In reality 
the chemical potential, and consequently the gap and the 
critical temperature, change with the distance from the 
center and the transition at the center and that at the 
edge of the core may not be simultaneously. For a typi- 
cal compact star, the density at the center is at most two 
times larger than that on the edge[l4j and the quark mat- 
ter cools very rapidly and reaches isothermal condition in 
a few hours because of the high thermal conductivity [l5| . 
As a crude estimation, we associate the chemical poten- 
tial of 500MeV to the center of the quark matter core and 
define its edge at the radius where the transition temper- 
ature drops by half. Then time scale for the transition at 
the center to that at the edge is about a few minutes to a 
few days depending on the direct Urea process [15]. This 
is much shorter than the typical age of a compact star 
(10 4 — 10 6 y), justifying our picture of the sudden release 
of the latent heat. The subject deserves further study 
to clarify more quantitatively the interplay between the 
phase transition and the structure and the evolution of a 
compact star. 
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